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We consider interacting ultracold fermions subject to Rashba spin-orbit coupling. We construct 
a single-branch interacting theory for the Fermi gas when the system is dilute enough so that the 
positive helicity branch is not occupied at all in the non-interacting ground state. We show that the 
theory is renormalizable in perturbation theory and therefore yields a model of polarized fermions 
that avoids a multi-channel treatment of the problem. Our results open the path towards a much 
more straightforward approach to the many-body physics of cold atoms subject to artificial vector 
potentials. 
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I. INTRODUCTION 

The generation of artificial gauge potentials for ultra- 
cold atoms by means of properly engineered laser fields 
[ll-Q constitutes a promising route towards the simula- 
tion of, for instance, Abelian 0] and non-Abelian ^ Q 
gauge theories. Interestingly, before we are able to reach 
that major goal, many other interesting systems have 
been proposed, such as the simulation Q of a bosonic in- 
teracting gauge theory @, while some others have even 
been already implemented experimentally. These include 
the generation of a uniform Abelian vector potential [lo| 
leading to the observation of quantized vortices in a Bose- 
Einstein condensate and spin-orbit coupling fl^ . 

Systems of ultracold spin-orbit coupled atoms have 
gained major interest recently. The phase diagram for 
weakly interacting bosons has been elucidated both in 
the untrapped case [Tsj , where the system exhibits plane- 
wave and striped phases, and in the case of a harmoni- 
cally trapped Bose gas [IJ], where the phase diagram is 
even richer. In the case of fermions, the strength of the 
spin-orbit coupling plays a major role in the formation 
of two-body bound states which exist even in the 
negative side of the s-wave scattering length in vacuum. 
Moreover, the crossover between a Bose-Einstein conden- 
sate (BEC) of molecules and a Bardeen-Cooper-SchriefFer 
(BCS) superfluid state can be driven by the spin-orbit 
coupling alone when the s-wave scattering length is fixed 
(for an extensive overview, see [l6j). 

The many-body problems with and without spin-orbit 
coupling are very different from each other. For spin- 
1/2 fermions, the spin-orbit term couples the two z- 
components of the spin, turning two- par ticle scattering 
into a genuine multi-channel process 17[. For a dilute 



Fermi gas, the non-interacting lowest-energy channel is 
strongly preferred. However, as interactions are turned 
on, a single-channel description breaks down as a non- 
trivial, seemingly non-renormalizable ultraviolet (UV) 
behavior appears [l8| . The full multi-channel descrip- 
tion of the many-body problem is free of such anomalous 
UV structure, as shown by Ozawa and Baym in [l7j . 



The multi-channel problem is a formidable one, hence, 
a renormalizable single-channel theory would greatly sim- 
plify the many-body problem. We here construct such a 
theory via perturbative renormalization of the two-body 
interaction, in a way that single-channel scattering re- 
produces the corresponding component of the exact T- 
matrix. We then obtain the energy of the Fermi gas 
in the normal phase to second order in the renormalized 
coupling constant, which is finite and independent of any 
momentum scales, thus showing that the theory is renor- 
malizable. 



II. SINGLE-PARTICLE SPECTRUM 

We begin by briefly reviewing the single-particle prob- 
lem, which also sets the notation for the rest of the article. 
The most general single-particle problem with spin-orbit 
coupling can be diagonalized exactly [l^. We restrict 
ourselves, for concreteness, to a particular system of in- 
terest, with spin-orbit coupling of the Rashba type. The 
single-particle Hamiltonian is given by 
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where a = {'yx,<^y,<^z) is the vector of spin-1/2 Pauli 
matrices. The helicity % = a ■ pj^, in the x — y plane, 
is defined as the component of the spin in the direction 
of the in-plane momentum pj_ = {px^py,0). The cor- 
responding eigenvectors \ip^^^) define what we call the 
helicity basis. Since [H, Hq] = 0, helicity is a good quan- 
tum number of the system, and we work, from now on, 
in the helicity basis. The two eigenvalues h of the helic- 
ity for momentum p^ are given by h± — ±p_\_, and the 
corresponding cigcnstates are given by 

where 7^ is the polar angle of pj^, given by tanjk ~ 
ky/kx- The eigenstates of the single-particle Hamiltonian 
arc therefore given by 

|V^(±)(r)) = |V^(±)(r^))e»'=^^ (3) 
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and the energy dispersions e± have the form 



e±(p) 



2m 
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. It is also convenient to define spin states of the form 

lt>±e'^Mi> 



± = 



V2 



(5) 



withfc^ = k^/kj_, so that = |k)(g)|±) = |k,±). In 

the following, we will refer to negative (positive) helicity 
fermions as lower-branch (upper-branch) fermions. 



III. FREE FERMI GAS 

Before we proceed with the inclusion of interactions, 
we study the ground state of the free Fermi gas. We are 
aiming at constructing a single branch theory. Therefore, 
we must be in the regime where the fermions in the non- 
interacting ground state only occupy the lower branch. 
Each fermion must therefore have an energy satisfying 



(6) 



where Ep is the Fermi energy. We define, for convenience 
Ep = kp/2m (with h = 1 throughout), although this 
definition is arbitrary, since the single-particle dispersion 
is not quadratic. The above condition gives the region of 
integration for the z-component of the momentum in 
the energy states 



(7) 



The integration interval for q±_ is given by {\—kp, \-\-kp)^ 
which is twice as large as the Fermi momentum. The 
q, = Q states therefore reach the energy Ep when q±^ 
equals either end of the interval. The ground state energy 
is then given by 
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dq±qi 
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After a tedious but straightforward integration, we ob- 
tain 



Ea 
V 



\k% 



where the density is given by 



47r 



(9) 
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It is interesting to note that the ground-state energy van- 
ishes, for large A, as Eq ^ 1/X in three-dimensions, while 
it would behave as 1/A^ in two dimensions. This is due 
to the fact that in 3D finite-density states need to acquire 
momentum in the z-direction, and therefore there are less 
states with energies close to zero as compared to the 2D 
case where virtually every occupied single-particle state 
has vanishing energy. 



IV. EFFECTIVE INTERACTION 

Several works on interacting bosons and fermions un- 
der spin-orbit coupling have noticed the appearance of an 
undesired ultraviolet (UV) logarithmic divergence in the 
T- matrix [l7| . both in the vacuum and in the medium, 
when all other helicity channels, which include intra- 
band interactions and interband interactions and tran- 
sitions, are not taken into account. As shown by Ozawa 
and Baym [l7| , the cancellation of the logarithmic diver- 
gence comes from the double transition (DT) process of 
two upper-branch fermions changing their helicity to (— ) 
(|-|--|-) — >■ I )). Excluding the DT process in the effec- 
tive theory means that the logarithmic divergence has to 
be explicitly eliminated. This can be done elegantly by 
introducing a pseudo- pot ential a la Fermi-Huang [20| in 
momentum space [2ll. |22|. All other collision processes, 
and the finite contribution of the DT process, are taken 
into account in the renormalization of the interaction. 



A. Logarithmic divergence 

The interaction between fermions in the lower-branch 
is given by 



2V ^ 

Pl+P2=P3+P4 



A (p) A (71, 72, 73, 74)c],,cJ,3Cp,Cpi , 

where Cp annihilates a fermion in state |p, — ), 5* is the 
bare (un-renormalized) coupling constant, and 



A(7i,72,73,74) 
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(12) 



The relative momentum p is defined by p = (p2 — Pi)/2, 
with 7i =_7pi, while A(p) is Tan's A-distribution |21 
given by |22| 



A(p) = 1 



<5(l/p) 



V 
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Note that by including Tan's distribution in the interac- 
tion in Eq. (fTTj) . the linear ultraviolet (UV) divergence 
connected to the s-wave scattering length is already ab- 
sent. The logarithmic divergence in the single-branch 
model comes from attempting to calculate the two-body 
T-matrix using the single-channel approach. To see this, 
we notice that the single-channel second Born approxi- 
mation is given by 

r*(p,p';Q) = 5,A(7i,72,73,74) 
■1 



9* 



dk A(k) ^(^1' '^2, 75, 76)A(75, 76, 73, 74) ^ ^^^3^ 



(2^) 



e_(f-k) + e_(§-fk) 



(14) 



In the above equation, Q is the center-of-mass momen- 
tum, and p and p' are, respectively, the incoming and 
outgoing relative momenta of the two-body system. The 



3 



angles 71, 72, 73, 74, 75 and 75 correspond, respectively 
to the momenta Q/2 + p, Q/2 - p, Q/2 - p', Q/2 + p', 
Q/2-k and Q/2+k . As it stands, Eq. ([14]) is free of the 
linear UV divergence (which is removed by the pseudopo- 
tential), but has a logarithmic divergence. After trivial 
algebraic manipulations, we obtain from Eq. (|14p 

T*(P,P';Q) = .9*A(7i, 72,73, 74)x 



and 



dk 

(2^ 



A(k)- 



2[1 - cos(75 
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The UV structure of the T-matrix at a large momentum 
cutoff /3, is given by 



5* A 



+ 
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where /3o is a finite, arbitrary momentum scale with the 
only purpose of rendering the argument of the logarithm 
dimensionless. The logarithmic divergence above is the 
one we aim at eliminating in a physically consistent way, 
that is, renormalizing it away without the appearance of 
any extra scales or fitting parameters in the system. We 
will show in the next section that this is indeed possible 
at the many-body level. 



B. Renormalization of the coupling constant 

The single-channel T-matrix, T* , calculated in the pre- 
vious subsection is still not correct at second order in 
the coupling constant. We therefore invoke perturbative 
renormalization so that the single-channel model repro- 
duces the correct T-matrix. This is given by p7j 



r(p,p';Q) = 



gA(7i,72,73,74) 



l + ^[A{Q/2) + B{Q/2)y 



(17) 
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The single-channel T-matrix in Eq. (jlSp is obviously cor- 
rect to first order in g. However, the second-order term 
is not correct and logarithmically divergent, and needs 
to be renormalized. In fact, the logarithmic divergence 
can be safely eliminated at this stage using a pseudopo- 
tential, as we will show in the last section of the article, 
but we first show that it is cancelled in the course of 
renormalization. We begin the renormalization process 
by expanding the bare coupling constant in powers of its 
renormalized counterpart 

,9*(A,g) = ,9 + aA,Q.9' + 0(5^). 
The renormalization condition reads 

T(p,p';Q)=r*(p,p';Q), 



(20) 



(21) 

to the given order in g. Expanding the exact T-matrix, 
Eq. (IT71) . to second order, and equating it to the single- 
channel T-matrix, Eq. (|15p . we obtain 
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where we have defined 
dk 



Ia,q — 2 



-A(k) 



-[A{Q/2) + BiQ/2)], (22) 
1 - cos(75 - 76) 



e„(Q/2-k)+e_(Q/2-|-k) 



(23) 

As can be readily checked, Ia.q, and therefore oa^q, is 
logarithmically divergent. 



where g = Ana/m is the renormalized coupling constant, 
with a the s-wave scattering length, Q ~ Q/ X and A and 
B are given by 



AiQ/2) 




V. INTERACTING FERMI GAS 

We develop here a perturbation theory for the inter- 
acting Fermi gas, using our pcrturbativcly-renormalizcd 
single-channel theory, and show that the theory is in- 
deed renormalizable and free of any logarithmic UV- 
divergence after renormalization. 

The first order contribution of the interaction is given 

by 




^ {F\V\F) 
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where \F) is the non- interacting ground state (Fermi sea) 
and the momenta in the integration are restricted to val- 
ues within the Fermi sea, together with g* = g + 0{g^). 
The final result is 



.9 2 



(25) 
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At second order there are two different contributions. 
The first one is due to the renormalization process at 
order in Eq. (|20l) . 



^_(9'\ 1 
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where Q = q -I- q' . The other contribution arises from 
the usual second-order diagrams and is given by 
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A(k)Ck,k'Cq^q' 



V A{2^fjp ^Jp ^ J^^p e_(k,k',q,q')' 

(27) 

where k' = q + q' — k, Ck,k' = 1 — cos(7k — 7k') and 

e„(k,k',q,q') = e_(k) +e_(k') - e_(q) - e_(q'). (28) 

There is still one important point yet to prove, and 
that is the independence of the results on the renormal- 
ization prescription, so that there is no arbitrary scale 
left. The theory will then be free of UV-divcrgcnccs and 
therefore has predictive power. We show here that the 
second order correction to the energy 
is finite by showing that the logarithmic divergences 



present in E\ ' and E)^' cancel each other. To see 
this, notice that for large k, we have l/e_(k, k', q, q') = 
l/2e_(k)-hO[(e_(k))-2]. All the divergent UV structure 
of the integral in Eq. ([27l) appears at this order alone. A 

(2) 

formal expansion for E2 then has the form 



p(2 



E. 



(2) 



V 



8(27r) 



.9 2 
—P 



A(k) 



(2) 



(29) 



while the divergent part of E\ ' is easily isolated from 
the finite part. This is given by 
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■p- / dk 



A(k) 

6_(k) 



(30) 



From the two equations above, we see that £''^' ~ 

(2) (2) 

E\ +E2 is not logarithmically divergent, which is what 
we wanted to show. The single-channel theory can there- 
fore be nicely renormalized in favor of only the scattering 
length. 



mic divergence and is finite at all steps. To do so, we 
need to replace the pseudopotential A(k) by a new pseu- 
dopotential A(k), together with an irrelevant, arbitrary 
momentum scale /3o- The modified pseudopotential has 
the form 



A(k)=A(k)- Aiog(^^)^(2)(i/fc^), 



(31) 



where 5'^'^^ is the two-dimensional Dirac delta. It is easy 



to see that in this case E^' and i?2^^ are both finite 
(2*^) separately, and so is the bare coupling constant p,, but 



both and E^' depend on the choice of /3o- The 

renormalizability of the theory depends now on proving 
that the results are independent of the scale /3o, that 
is, to show that it is irrelevant. Proceeding in a way 
completely analogous to that of the previous section, we 
find that this is indeed the case. 
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VII. CONCLUSIONS AND OUTLOOK 



We have introduced a renormalizable theory for in- 
teracting fermions subject to spin-orbit coupling. Our 
theory is valid in the dilute regime where the fermions 
in the non-interacting ground state occupy only the 
lower hclicity branch. The effective single-branch model 
corresponds to interacting polarized fermions, and thus 
opens the path to a simpler treatment of the many-body 
problem, circumventing the intricacies of the full multi- 
channel system. We have illustrated our methods by cal- 
culating the second-order correction to the ground-state 
energy of the repulsive Fermi gas. 

As a natural extension of our work, it would be inter- 
esting to calculate the non-Hermitian optical potential, 
which takes explicitly into account the population of the 
other channels. Our model may also be of interest for a 
more sofisticated treatment of the BEC-BCS crossover, 
and for reduced dimensional Fermi and Bose gases sub- 
ject to spin-orbit coupling. 
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